Rabi-type oscillations in a classical Josephson junction 
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We study analytically and numerically the phase-modulation properties of a biased classical 
Josephson tunnel junction in the zero- voltage state and phase-locked to an external ac field. We show 
that the phase-locked state is being modulated in the transients, or in response to perturbations, 
and the modulation frequency is calculated as a function of relevant system parameters, such as 
microwave field amplitude. The numerical analysis is parameterized similarly to recent experimen- 
tal results in which a combination of pulsed ac signals and relative switching from the zero-voltage 
state are used to probe the internal excitations of the junctions. Our analysis demonstrates that 
the modulation of a phase-locked state in an entirely classical Josephson junction produces oscilla- 
tions analogous to quantum mechanical Rabi oscillations, expected to be observed under the same 
conditions. 



The symmetric oscillations of a quantum system be- 
tween its ground and excited states have recently at- 
tracted much interest within the field of quantum coher- 
ence and quantum computation; see, e.g.ji. Existence of 
such oscillations was originally reported by Rab£, who 
studied the time dependent dynamics of a spin subject 
to an external time-varying magnetic field. The phe- 
nomenon is presently considered as a tool in characteriz- 
ing coherence (and decoherence) of quantum states. Due 
to the possibility offered by superconducting circuit inte- 
gration, recent attention has been devoted to detect Rabi 
oscillations in macroscopic quantum systems based on 
Josephson junctions. The reported experiments^ indi- 
cate evidence of phase oscillations in Josephson junctions 
similar to what can be expected from the transitioning 
between the quantum states of a multi- level atom. In this 
letter we demonstrate that the energetic properties of dif- 
ferent dynamical states in a nonlinear classical Josephson 
junction can produce phase oscillations with observable 
features analogous to the reported quantum phenomena. 

The above mentioned investigations of Rabi 
oscillations^ were concerned with excitations of 
the zero-voltage of Josephson tunnel junctions driven 
by weak external microwave signals. The properties 
of dc and ac driven Josephson junctions have been 
investigated in several contexts^, but most of the at- 
tention has been devoted to the phase-locking of nonzero 
dc-voltage states such as constant-voltage Shapiro steps 
in the current-voltage characteristics. Our investigation 
is entirely focused on the phase-locking of a Josephson 
junction in the zero-voltage state. We assume that the 
relevant equation of motion for this system is given 
by the classical RSCJ (Resistively and Capacitively 
Shunted Junction) models, 



(p + aip + sin ip = r) + e s sin oj s t . 



(1) 



Here sini^, r\, and £ s sintj s < represent Josephson, dc- 
bias, and ac-bias currents, respectively, normalized to the 
Josephson critical current I c . The amplitude and temper- 



ature dependence of I c is established by the Ambegaokar- 
Baratoff equation^. Time is normalized to the in- 
verse of the Josephson plasma frequency Wq 1 , where 
u% = 2eI c /hC = 2irI c /$ C, $ = h/2e = 2.07-10~ 15 Wb, 
and C the capacitance of the junction. Damping is rep- 
resented through a = fnu /2eRI c — $o^o/27r RI c , where 
R is the model resistance. Normalized energy is defined 
by 



2^ 



1 — COS ip — r)(f 



(2) 



where the characteristic energy is Hj = I c <&o/2tt, and 
the resulting input power is 



H = e s ipsmoj s t — aip 



(3) 



Our analysis of phase-locking in the system described 
above take into account the anharmonicity of the po- 
tential through the zero- voltage, phase-locked monochro- 
matic ansatz, 



ip = ipo + tp = ipo + a sm(u> s t + t 



(4) 



where (ip) = (po = 0, and a is an oscillation amplitude. 
Inserting this ansatz into Eq. (1) gives the static, 



Jo (a) simpo = t] 



(5) 



component, where J n is the Bessel function of the first 
kind and nth order. The dynamic component is deter- 
mined by 



t/j + aip + cos po sin ip = e s sin uj s t 



(6) 



Inserting (4) and (5) into (6) gives the steady-state non- 
linear relationships, 
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FIG. 1: Comparison between numerical results (•) and ana- 
lytical ansatz Eq. (10) (solid lines) for anharmonic relation- 
ships between (a) frequency and amplitude Eq. (9), (b) aver- 
age phase and amplitude Eq. (5), and (c) energy and ampli- 
tude Eq. (11). Data shown for rj = 0.94259 and a = e a = 0. 



where uo r is the anharmonic resonance frequency 
2 _ 2Ji(a) 
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(9) 



Notice that uj r is the natural anharmonic frequency of an 
excitation with amplitude a when e s = a = 0. Thus, the 
solution (for e s = a = 0) establishes a useful ansatz, 



tp — ipo + a sm(uj r t + i 



(10) 



for the unperturbed problem if tpa is given by Eq. (5). 

We substantiate this very simple ansatz through Fig. 1, 
which compares direct numerical simulations of Eq. (1) 
with the predicted relationships (9) and (5), as well as the 
energy H — Hq, where Hq = 1 — \f\ — rj 2 — rj sin -1 77, as 
a function of a. The average total energy for the ansatz 
Eq. (10) is 
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(11) 



Specific system parameters are: r\ = 0.94259 and a = 
e s = 0. We clearly observe that the ansatz represents 
the anharmonicity very well for small to moderate am- 
plitudes. 

Figure 2 shows the comparisons between the pertur- 
bation analysis and numerical simulations for the phase- 
locked states when a = 10" 



T) = 0.94259 for two differ- 
2 - 0.577886 (Fig. 2a,b) and 



cnt frequencies uj s = \J\ — r\ 
lo s = 0.5237246 (Fig. 2c, d), the latter being the resonance 
frequency for an amplitude close to the switching point 
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FIG. 2: Oscillation amplitude a and energy (H) — Ho as a 
function of microwave a mplitude e s for a = 10 -3 and signal 
frequencies (a-b) u) s = ij\ — rj 2 — 0.577886, and (c-d) u s = 
0.5237246. Numerical simulations of Eq. (1) are represented 
by markers (•), and solid lines represent Eqs. (7) and (12). 



from the zero-voltage stated. We clearly see the close 
agreement between theory and simulations outlining the 
anharmonic saturation of the resonance in Fig. 2a,b, and 
the switching between small amplitude (off-resonance) 
and large amplitude (resonant) excitations in Fig. 2c, d. 
The average energy for the phase-locked system is 



H = -a 2 ^ + 1 
4 



,aJ (a) 
2Ma) 



77 sin 



-if V 



\Jo(a) 



(12) 



Thus, the phase-locked states seem adequately described 
by the anharmonic perturbation method outlined above. 
We will later use uj s — 0.577886 as the signal and oj p — 
0.5237246 as the probe frequency for demonstrating the 
classical equivalent to Rabi oscillations. 

Stability of the phase-locked state can be evaluated by 
following the treatment outlined for breathers in driven 
sine-Gordon systems^. The energy change during one 
microwave period is 



AH 



Hdt 



(13) 



-e s ansm8 - ana 2 u s = I {ln) sin - I {out) , (14) 

where Eq. (14) is the result for the monochromatic 
ansatz, and where steady state phase-locking is charac- 
terized by, 



AH 







sin Uq 



j(out) 
I (in) 



(15) 



Notice that this power-balance expression is consistent 
with (7) and (8). 

A perturbation to a phase-locked state may introduce 
slow (transient or steady state) modulation frequencies to 
the dynamics. Such modulations have also been observed 
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FIG. 3: Frequency/amplitude and energy modulation as a 
function of time after the onset of a microwave field with 
amplitude e s = 10 -3 . Other system parameters are a = 10 -3 , 
uj s = \/l — rj 2 , and rj = 0.94259. Results are obtained from 
Eq. (18). 



and described in the context of breather stabilization in 
perturbed sine-Gordon equations 8 . We give specific re- 
sults below for the phase-locked RSCJ model. 

We assume that a modulation frequency to the phase- 
locked state is slow and that the amplitude is small, i.e., 
that the modulated phase is given by 8 = 8q + 86, with 
1 50 1 <C 1 and \S0\ <C w«. In this approximation, we can 
write the change in total energy over one driving period 
as 

AH = = = iWcobOoM 

uo s uj s oa ouo r 

2 8H da ■■ 



lo s da doj r 



a^e 2 s - (aau s ) 2 60 . (16) 



Equation (16) provides a slow modulation frequency Q,r 
given by 



1 



dujr (dH\ 



ril = -- 2 . s aVe>-(aa. s y^^) (17) 

where a and £ s are related through Eq. (7), and a is 
related to uj r through Eq. (9). Notice that doj r /da < 0, 
leaving Eq. (17) with a (real) positive value. 

This modulation can be directly observed as transient 
behavior for t > in the equation 



Cp + aip + sin ip — rj + Q(t)e s sin(u; s i + S ) 



(18) 



where Q(t) is Heaviside's step function, and S is a con- 
stant, but random, phase. Figure 3 shows an example for 
a = 1(T 3 , T] = 0.94259, u s = 0.577886, and e s = 10~ 3 . 
We observe clear evidence of the slow modulation fre- 
quency £Ir around the phase-locked state at oj s of phase, 
amplitude, and energy. This behavior is present regard- 
less of the value of S . Thus, for any value of 9 S we can 
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FIG. 4: Dynamic Modulation as a function of microwave 



amplitu de e 3 . Parameters are: rj = 0.94259, a 



10" 



= {/l - rj 2 = 0.577886. (a) Modulation frequency Q R . 
Numerical simulations of Eq. (18) are represented by markers 
(•) and the solid line represents Eq. (17). (b) Switching prob- 
ability. Each point represents 10 4 events simulated through 
Eq. (20) with parameters: uj v = 0.5237246, e v = 0.004, 
k B T/Ej = 5 ■ 10~ 5 , tl = 10 4 , ti=tl + 4043, tj = t\ + 3185, 
tp = tl + 6248. All switching events take place in the inter- 
val ]tl;ti]. (c) Variation in switching probability. Solid line 
represents the variation as derived from Eq. (17). Markers 
represent simulations of Eq. (20), analysed through Eq. (21), 
with parameters and microwave timings tj — tl and tp — tt set 
as in (b), but with probe onset times given by At 3p — 2205 
(A), At sp = 3185 (o), At sp = 4165 (x), and At sp = 5145 
(□). 



expect that the energy is at its minimum Hq for t < 0, 
and oscillates around H for t > with frequency f2#. 
Combining this information, we can now write the time 
evolution of the energy as 



Hit) 



H 

H-(H- H )e 



, t < 
^'cosfiflt , t > 



(19) 



where /3 1 is the transient timeifl. 

Figure 4a shows a direct comparison between numeri- 
cal simulations and the perturbation result Eq. (17) of the 
modulation frequency as a function of the microwave am- 
plitude £ s . The overall agreement is very good, and the 
discrepancy is understandable considering the monochro- 
matic ansatz, which neglects energy contributions from 
higher harmonics. 

The above analysis can be applied directly to recent 
experimental reports 3,4 of Rabi oscillations. A complete 
system with signal and probe is governed by 



(p + aip - 

+mi) 



simp = -q + (6(4) 
- Q(t*))e p sin(w p i + 



- Q(t l s ))e s sm(u! s t - 
%) + n{t) , 



(20) 
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where t\ < t\, t\ < t p < and n(t) is a ther- 
mal noise term determined by the fluctuation-dissipation 
relationship 11 . Following the experimental procedure, we 
choose the probe frequency smaller than the signal fre- 
quency, ui p < oj s w y/l — rj 2 , and e p is chosen such that 
phase-locking at ui p presents two possible energy states 
(see Fig. 2c, d), one small amplitude and one near or be- 
yond switching. Switching is then induced with high 
probability if the probe field is causing a high energy 
state at u) p , while low probability of switching occurs if 
the probe field excites the low energy state (see Fig. 2d). 
The determining factor is the timing At sp = t 1 — t 1 be- 
tween the onset of the two microwave fields. This is ap- 
parent from Fig. 3 and Eq. (21) since UfjAt sp — 2np (p 
being an integer) will provide for a relatively low onset 
energy for the probe field with subsequent low probabil- 
ity for exciting a high-energy state at u> p . However, for 
fliiAt sp — 2np + 7T, the onset energy of the system is 
relatively large, and a much higher probability for high 
energy excitation at lo p can be expected. A reasonable set 
of probe parameters can be chosen to ui p — 0.5237246 and 
e p = 0.004. As mentioned above, this probe frequency 
provides a resonant amplitude that takes the junction 
phase close to the saddle point of the energy potential, 
and the amplitude e p is chosen such that both high and 
low energy w p -states can be chosen based on the state of 
the system at time t = tt. This is the recipe for the clas- 
sical equivalent of Rabi oscillations in microwave driven 
Josephson junctions, and they can be observed in several 
different ways. One is to keep e a constant while varying 
At sp for the direct observation of the modulation fre- 
quency as has been done inl Another is to keep At sp 
constant while varying e s , as was demonstrated in 3 . We 
will demonstrate the latter in Fig. 4b, where we show sim- 
ulation results for At sp « 3185. Each point on the figure 
represents 10,000 events of sequentially applying the sig- 
nal and the probe, each with random phases S and 9 P , 
and detecting whether or not the system switched to the 
non-zero voltage state within the probe duration t p — t p . 
The fraction of switching events is shown in Fig. 4b. We 
note that this set of data has been generated with a small 
amount of thermal noise (ksT/Ej w 5 • 10 -5 ) in order to 
mimic experimental reality. However, this is not signif- 
icantly affecting the result, since the random phases 6 S 
and 9 p are the primary factors in determining if a switch- 



ing event occurs. The variations in switching probability 
as a function of e s are very clear from this figure, and the 
characteristics of the observed variation is given by 



where T is defined in Fig. 4b. A convenient way of quanti- 
tatively comparing the modulation in the switching prob- 
ability as a function of microwave amplitude e s to the 
predicted behavior related to the modulation frequency 
£Ir is to show r/ At sp = dVtn/de s as a function of e s for 
different At sp . This is done in Fig. 4c, where the solid 
curve is derived from Eq. (18), and the markers represent 
simulation results as shown in Fig. 4b. Good consistent 
agreement is found between the classical theory and the 
modulation of the probability of switching in the range 
where we have obtained data. The trend in the compar- 
ison, i.e., all the simulation results being above the the- 
ory line, is directly related to the monochromatic ansatz, 
which underestimates the modulation frequency (Fig. 4a) 
and the system energy (Fig. 2) of the phase-locked state. 
The above analysis and simulations demonstrate a di- 
rect analogy between transient modulation in a classical 
nonlinear Josephson junction irradiated with microwaves 
and reported Rabi oscillations in a system under similar 
conditions. 

In conclusion, we have shown that energy level degen- 
eration in the nonlinear system describing a Josephson 
junction can lead to transient oscillations when the phase 
of the junction is modulated by adequately shaped exter- 
nal ac-current terms. The resulting dynamical features 
demonstrate close analogy to the quantum phenomenon 
of Rabi oscillations. We note that a similar analogy be- 
tween classical and quantum analysis in relation to Rabi 
oscillations and splitting has been investigated in quan- 
tum optics^. As far as Josephson effect is concerned, 
analogies between classical and quantum descriptions of 
experimental results have been demonstrated by recent 
comparisons of theoretical and experimental results*^ 
relating to resonant switching from zero- voltage states. 
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